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Abstract. We study two models consisting of reflecting one-dimensional Brownian 
"particles" of positive radius. We show that the stationary empirical distributions 
for the particle systems do not converge to the harmonic function for the generator 
of the individual particle process, unlike in the case when the particles are infinitely 
small. 



1. Introduction 

In this article we consider the dynamics of a collection of hard Brownian spheres 
with drifts or boundary conditions that includes instantaneous reflections upon col- 
lisions. The models are similar to existing ones in the literature that consider point 
masses instead of spheres of a positive radius. We will show that the (empirical) 
distribution of a family of Brownian spheres behaves differently from the (empirical) 
distribution of the point Brownian particles in some natural models. In particular, 
the distribution of Brownian spheres fails to satisfy the usual heat equation under 
circumstances that lead to the heat equation for the infinitely many infinitcsimally 
small Brownian particles. 

Various models of colliding Brownian particles have been considered in the sta- 
tistical physics literature. One stream, pioneered by Harris (1965), considers a 
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countable collection of Brownian point masses on the line that collide and reflect 
instantaneously. Also see the follow-up work on tagged particle in the Harris model 
by Diirr, Goldstein, & Lebowitz, Durr et al. (1985). A variation on the theme has 
been to replace the instantaneous reflection by a potential, and goes by the name 
of gradient systems. In these gradient systems, one studies the behavior of count- 
ably many particles under a repelling potential. Usually the potential is modeled 
as smooth with a singularity at zero; see the article by Cepa and Lepingle (1997). 
A particular example of this class includes the famous Dyson Brownian motion 
from Random Matrix theory; see Dyson (1962), and Cepa and Lepingle (2001). 
The other class of models, closer to our article, goes by the name of hard-core in- 
teractions, in which the Brownian particles are assumed to be hard balls of small 
radius, and consequently, there is instantaneous reflection whenever two such balls 
collide (plus possible additional interactions) . This is the spirit taken in the articles 
by Dobrushin & Fritz in dimension one, Dobrushin and Fritz (1977), and Fritz & 
Dobrushin in dimension two, Fritz and Dobrushin (1977), Lang (1977a,b) (with a 
correction by Shiga, 1979). The main focus of these authors is the non-equilibrium 
dynamics of the gradient systems. Also see the articles by Osada (1996, 1998), and 
Tanemura (1996) all of which consider properties of a tagged particle in the infinite 
system. 

In the discrete case, the various models of symmetric and asymmetric exclusion 
processes have been considered. Closest in spirit to the models discussed here is 
the totally asymmetric exclusion process (TASEP) considered by Baik, Dcift, & 
Johansson, Baik et al. (1999), and Johansson (2000) in connection with random 
matrices and the longest increasing subsequence problem. Specifically, if the initial 
configuration in TASEP is Z~, then the probability that a particle initially at —m 
moves at least n steps to the right in time t equals the probability distribution of 
the largest eigenvalue in a unitary Lagucrre random matrix ensemble. In recent 
subsequent articles Tracy and Widom (2008, 2009), Tracy & Widom explicitly 
compute transition probabilities of individual particles in the asymmetric exclusion 
process, extending Johansson's work. 

In this paper, we consider only one dimensional models, so our "spheres" are 
actually intervals. The title of this paper reflects our intention to study multidi- 
mensional models in future articles. We leave more detailed discussion to Section 4. 
That section also contains references to related research projects. 

We consider two models, which have the following common features. Informally 
speaking, both models consist of families of Brownian "particles" . The fc-th "par- 
ticle" is represented by an interval I k = (X k — e/2,X t k + e/2), where X t k is a 
Brownian-like process. The intervals I k and P arc always disjoint, for k ^ j. The 
processes X k are driven by independent Brownian motions. When two intervals 
I k and P collide, they reflect instantaneously. In the first model, the number of 
particles is constant and they are pushed by a barrier moving at a constant speed. 
In the second model, particles enter the interval [0, 1] at the left, they reflect at 
0, and they are killed when they hit the right endpoint. The second model is our 
primary focus because it is related to other models considered in mathematical 
physics literature — see Section 4. 

We are grateful to Thierry Bodincau, Pablo Ferrari, Claudio Landim, Mario 
Primicerio and Jeremy Quastcl for very helpful advice. We would like to thank the 
referee for the suggestions for improvement. 
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2. Extreme crowding 

We start with an informal description of our first model, which consists of a 
fixed number n of "particles" . The fc-th leftmost "particle" is represented by an 
interval/* = (X k — e/2, X k + e/2). The intervals I k and P are always disjoint. The 
processes X k are driven by independent Brownian motions. When two intervals I k 
and P collide, they reflect instantaneously. The intervals are pushed from the left 
by a barrier with a constant velocity, that is, the leftmost interval reflects on the 
line x = ct. 

Formally, we define {X°, X 1 , . . . , X n } to be continuous processes such that X® = 
-e/2 + ct, X k - X^ 1 > e for all k > 1 and all t > 0, and for k > 1, 

dX k = dB k + dL k - dM k , 

where {B 1 , . . . , B n } arc iid Brownian motions, and L k and M k are nondecreasing 
processes such that 

/>oo />oo 
J 1 {X*-X*- 1 >e} dL t = and J 1 {X^ 1 -X*>s} dM t = 0. 

(Here, we may interpret X n+1 = oo.) The distributions of Xq for 1 < k < n will 
be specified later. 

To construct the solution to this Skorohod problem, consider first the processes 
Y k = X k -(k- l)e - e/2 - ct. These processes satisfy Y° = 0, Y k - Y^ 1 > for 
all k > 1 and all t > 0, and for k > 1, dY k = dB k -cdt + dL k - dM k , where 

/>oo 

l{Y t *-Y' t ' ! ~ 1 >o} dL t = and y ^{Y t k+1 -Y t k >o} = 0- 

We may therefore construct the processes {Y 1 ,...,Y n } using order statistics. 
Namely, let {Z 1 , . . . , Z n } be defined by dZ k = dB k - cdt, and reflected at 0. For 
every fixed t > 0, we let F/, Y t 2 , . . . , Y t n be ordered Z t fc 's, that is, {F/, . . . , F™} = 
{Z t x , . . .,Z?} andFj 1 < F t 2 < • • • < F™. Finally, we let X k = Y k + (k-l)e+e /2 + ct 
&ndI t k = (X?- S /2,X t k + E/2). 

Let ne = b. We will fix b > and analyze the behavior of the system of intervals 
{I k } as n — > oo. In other words, we will keep the total length of all intervals I k 
constant. 

The stationary distribution for Z k has the density <p{z) = ce~ cz for z > 0, with 
c = 2ci, because 

1 d 2 . . d , . 
2^ (Z) + C1 ^ )= °- 
Consider any < xi < x 2 < oo, let A denote the Lebesgue measure, and let 

x({xi+ c t lX2 + ct]n\J 1<k<n i k ) 

d([xi,x 2 ]) = d t ([x 1 ,x 2 ]) = — ^ — '-. (2.1) 

x 2 - x\ 

The quantity d([xi, Xi\) represents the average density of "particles" I k on the 
interval [x\, x%\. 

We will say that the intervals {/ fc } have the pseudo-stationary distribution if 
all Z k, s arc independent and have the stationary distribution ip for t = and, 
therefore, for every t > 0. 
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Theorem 2.1. Suppose that the intervals {I k } have the pseudo-stationary distri- 
bution. Fix arbitrary p\,d\ < 1, d<z > 0, and < X\ < x 2 < b < X3 < X4 < 00. 
There exist Co, no < 00 such that for c > cq, n > uq and any t > 0, we have 

P{d t ([x 1 ,x 2 })>d 1 )>p 1 , (2.2) 
P{d t {[x 3 ,x 4 }) <d 2 )> Pl . (2.3) 

The theorem says that the "particles" I k clump together and there is a sharp 
transition in density of "mass" around x — b. This is in contrast with infinitely 
small "particles" Z k whose empirical distribution is close to the distribution with 
the density ip(z) = ce~ cz that displays no sharp drop-off. 

Proof of Theorem 2.1: Without loss of generality, we let t = 0. We define y k £ 
(0, 00) in an implicit way by the following formula, for k = 1,2, 3, 4, 

rvk 

x k = b ip(z)dz + y k . 
Jo 

Note that y\ < y 2 , and that for e > sufficiently small (that is, for n = 6e _1 
sufficiently large), it is possible to choose 2/5,2/6 such that y\ < 2/5 < 2/6 < 2/2, and 

b f VB (p(z)dz - 2s b f y ' 2 V(z)dz 

2/2-2/1+0 J yi <p{z)dz 2/2-2/1 + J yi ip(z)dz 

Since 6 — x 2 > 0, we can find c so large that, 

c(b - x 2 ) 



1 + c(6 - x 2 ) 



>l-(l-di)/2. (2.5) 



Let [a] denote the smallest integer greater than or equal to a. By the law of large 
numbers, if n is sufficiently large, the number of Zq 's in the interval [0, 2/1] is smaller 
than or equal to n (p(z)dz, with probability greater than 1 — (1 —pi)/2. If this 
event holds then there are exactly |~n J Q y5 <p{z)dz~\ processes Zq in some (random) 
interval [0,2/7] with 2/7 > 2/1 • This implies that there are exactly \n J^ 5 <p(z)dz~\ 
processes Xq in [0,e \n J^ 5 ip(z)dz~\ + y 7 ]. Note that 



2/5 

n I (p(z)dz 




rV5 ryi 

2/7 > 6 / (p(z)dz + 2/7 > b / Lp{z)dz + yi = xi. 



Hence, the number of Xg's in the interval [0,xi] is smaller than or equal to 
n Jq 5 (p(z)dz + 1, with probability greater than 1 — (1 — pi)/2. A completely anal- 
ogous argument shows that, if n is sufficiently large, then the number of Xq 's in 
the interval [22,00] is smaller than or equal to n J ip(z)dz + 1, with probabil- 
ity greater than 1 — (1 — pi)/2. Both events hold with probability greater than 
1 — 2(1 —px)/2 = pi, and then the number of Xq's m [^l, 2 ^] is greater than or 
equal to n ip(z)dz — 2. This and (2.4) imply that 

en tp{z)dz -2s b tp(z)dz - 2e 

rV2 



d{[x 1 ,X 2 \) > 



X2-X1 b J Q V2 ip{z)dz + 2/2-6 J Q yi ip(z)dz - 2/1 



b f V6 ip(z)dz - 2e b f V2 <p(z)dz 

V5 m , ^ > m Pm , ; , - (1 - di)/2. (2.6) 



2/2-2/1 + 6 III (p(z)dz 2/2-2/1+6 III <p(z)dz 
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We have 

x 2 = b I (p(z)dz + y 2 = b I ce~ cz dz + y 2 = y% + b - be~ cy2 , 



so e cy2 = (y 2 ~ x 2 + b)/b and, therefore, for z < y 2 , 

ip{z) = ce~ cz > ce~ CV2 = (c/b)(y 2 - x 2 + b). 

We combine this estimate with (2.6) and (2.5) to see that, with probability greater 
than pi , 



b f V2 (p(z)dz 

d{[x u x 2 ]) > \ ' ' (1 - dO/2 

V2 - Vi +oj yi <p(z)dz 



> tJZWMy* -*> + <>)** _ (1 _ di)/2 

V2 - yi +bj^(c/b)(y 2 -x 2 + b)dz 

= c (^2 -Vi)(V2 -x 2 + b) _ (l - dx) 12 
V2-y\+ c(y 2 - Ui){y2 -x 2 + b) 
c{y 2 -x 2 + b) 



1 + c{y 2 -x 2 +b) 



1 + c(o - araj 
> 1 - (1 - dx)/2 - (1 - dx)/2 = d x . 

This completes the proof of (2.2). The proof of (2.3) is completely analogous. □ 



3. Brownian gas under pressure 

In this model, "particles" I k arc confined to the interval [0,1]. More precisely, 
their centers are confined to this interval. The k-th leftmost "particle" is represented 
by an interval I t fe = (X^—e/2, X^+e/2). The intervals I k and P are always disjoint. 
The processes X k are driven by independent Brownian motions with the diffusion 
coefficient a 2 . When two intervals I k and P collide, they reflect instantaneously. 
The particles are added to the system at the left endpoint of [0, 1] at a constant 
rate. In other words, they are pushed in at the speed a, so that a new particle enters 
the interval every e/a units of time. As soon as X k reaches 0, it starts moving as 
a Brownian motion reflected at 0. The fc-th interval is removed from the system 
when X k hits the right endpoint of [0, 1]. 

Formally, wc define {X 1 , X 2 , . . .} to be a collection of right-continuous, [0, oo]- 
valued processes such that 

X fe = -ke + e/2 for all k, (3.1) 
If S k = inf{i > : X k _ = 1 - e/2}, 

then X k is continuous on [0, Sk) and X k = oo for all £ > Sk, (3-2) 

X k _ x k+i > £ for all fc > x and all t > 0; and (3.3) 

dX k = { adt if £ £ [0,fce/a), 

* 1 a dB k + dL k - dM k if £ G [ke/a, S k ), 
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where a and a are positive constants, {B 1 ,!! 2 , . . .} are iid Brownian motions, and 
L k and M k are nondecreasing processes such that 

1 {X*-X?+ 1 >e} dL t =° and / 1 {X*- 1 -X(>e} dM t =0 - 
ke/a J he I a 

(Here, we may interpret X° = oo.) 

To construct the solution to this Skorohod problem, consider first the processes 
Y k = X k + ke — e/2 — at. These processes satisfy 

Y k = for all k, (3.5) 
If S k = inf{< >0:Y t k _ = l-at+ {k - l)e}, 

then Y k is continuous on [0, Sk) and Y k = oo for all t > Sk, (3-6) 

yk _ Y k+i > q for a ii > i an d a u £ > 0, and (3.7) 

dY k = l° ift€[0,ks/a), 

4 I a dB k -adt + dL k - dM k if t e [he /a, S k ), ' 



where 



Sk pS, 

rk 



lr-^fc v t+is.„i rfif = and / l, v t-i vk^^dM, 



{ yk_yk + l >0} U^ t -V O^U I ^{ Y -*- 1 -y t »>0} 



ke / a J ke I a 



Again, we shall construct the processes {F 1 ,^ 2 , . . .} using order statistics. 

Let Z k be a [0, oo)-valued process, satisfying the SDE dZ k = adB k — adt, and 
reflected at 0. The process Z k is defined on the time interval t S [fcsr/a, oo), and 
starts at Z k s / a = 0. At anytime i £ [ke/a, (k+l)e/a), only processes Z^ , 1 < j < fc, 
are defined. Let [aj denote the greatest integer less than or equal to a, and 

50 = 0, 

•4j = {ieZ:l<j<[ta/eJ}. t>0, 

51 = inf{< > : sup Z\ > 1 - at}, 

A h t = At 1 \ {m G Z : Zf k _ t = 1 - at + (k - 2)e}, t > S k -i,k > 2, 
5 fe = inf{i > 5 fc _i : sup Z{ >l—at+(k — l)e}, k > 2. 

Note that it is possible that A k = for some random k and t > 0. 
Convention ( C). For the sake of future reference, it is convenient to say that the 
process Z m is killed at the time Sk-i, where {m} = Ag ^ \ -4| fc i . In other 
words, Z m is killed when the corresponding interval / , defined below, hits the 
right endpoint of the interval [0, 1]. 

For every t G [Sfc-i, Sk), note that there are [ta/e\ — (k— 1) elements in A\. Let 
Y k ,Y k+1 , Y t lta/el be reverse-ordered Z{'s, j e A\, that is, {Y k , . . . , F t Lta/eJ } = 
{Z t J , j e A fc } and > ^t fe+1 > • > F t Lta/£j . Let Y t j = oo for j < k and F/' = 
for j > [ta/e\ . 
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It is elementary to check that {Y 1 , Y 2 , . . .} satisfy (3.5)- (3.8). We may therefore 
define 

X^ = Y t k - ke + e/2 + at, 
I k = (X k -e/2,X k + e/2). 

We have to modify slightly the definition (2.1) of density to match the current 
model. For t e [Sk-i, Sk), let 

A ([xi,x 2 ] nUfc<i<ito/ s j 7 * ) 

d {[x u x 2 ]) = d t ([x 1} x a =— 1 ■ (3-9) 

x 2 - Xl 

Theorem 3.1. Fix arbitrary < x\ < x 2 < 1, Pi < 1 and a, a, cq > 0. There exist 
to < oo and Eq > smc/i i/iai /or t > to and e G (0, eo), 

P fl 1 ~ X t /(r) , -c <d t ([x 1 ,x 2 ])<- 1 —^——+ Co ) > Pl . (3.10) 

\1 — x 2 + a*/\2a) 1 — x\ + a*/ [2a) J 

Intuitively speaking, the theorem says that the mass density at x € (0, 1) is close 
to (1 — x)/(l — x + a 2 /(2a)), for large t and small e. 

Proof of Theorem 3.1: We will use the coupling technique. Recall processes 
Z 1 ,Z 2 , . . . used in the definition of y fc 's — we will use the same Z k, s to construct 
auxiliary processes. Fix some v\ > 0, let Sk = infjt > : Z\ = Vi}, and let Z t fe be 
the process Z k killed at the time Sk ■ Let nt be the number of processes Z k alive at 
time t. Let Y^Y 2 ,..., Y t nt be ordered Z^s, that is, {Y?, . . . , Y t nt } = {Z 3 t ,S 3 > t} 
and < Y t 2 < ■ ■ ■ < Y t nt . For t e [S k -i,S k ) and j = k, . . . , k + n t - 1, wc let 

XI = Y? t+k ~i + (n t + k-j- l)e + e/2 + (t - [ta/e\e/a)a, 
Ti = (Xi-e/2,Xf+e/2). 

Every process X\ is defined on the interval [je/a,Sj) and it is continuous on this 
interval. Although it may not be apparent from the above formulas, the processes 
F J , X^ and P are constructed from 's in the same way as Y^ , X^ and P were 
constructed from Z^'s. We leave the verification of this claim to the reader. 

We will find the Green function G Vl (v) of Z k , i.e., the density of its occupation 
measure. Consider a process V with values in [— vi,Vi], satisfying the SDE dVt = 
dB t — a sign (Vt) dt, where B is Brownian motion, Vq = 0, and such that V is killed 
when it hits —v\ or v\. Note that the Green function G] ( «') of V is oik- half of 
G Vl (v) for v > 0. The scale function S(v) and the speed measure m(v) for V can 
be calculated as follows (see Karlin and Taylor, 1981, pp. 194-195), 



s(v) 


=exp U { 






S(v) 


Pt) 

= / s(x)dx = 




Jo 


m(v) 





2asign(x)/cr 2 )(ia; ) = cxp(2ausign(u)/cr' :! ), 
sign(w)cr 2 



2a 



(cxp(2aw sign(u) / a 



We will use formula (3.11) on page 197 of Karlin and Taylor (1981). In that formula, 
we take x — 0, so u(0) = 1/2, by symmetry. We apply the formula to functions g(v) 
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of the form g(v) = lr J , 3j „ 4 ] (w), to conclude that for v G (0,i>i), the Green function 
G v ± ( v ) is gi ven b y 

Gl(v) = (S( Vl ) - S(v))m(v) 
1 



— (exp(2awi/cr 2 ) — exp(2av/cr 2 )) exp(— 2av/a 2 ) 
2a 

1 



It follows that 



2a (exp(2a(w 1 -t;)/ ( T 2 )-l). 



G Vl (v) = 2Gj» = -(exp(2a(u 1 - v)/a 2 ) 
1 a 



Define vo G (0, oo) by setting 

ip(v) = a( 
and the following condition 



(p(v) = aG Vo {v) = i(exp(2a(w - v)/cr 2 ) - 1), 



wo 

2\ 



1= tp{v)dv + v = / (cxp(2a(u - v)/a 2 ) - l)dv + v (3.11) 
Jo Jo 

= {a 2 /2a)(exp(2av /(T 2 ) _ i). 
We define y k G (0, oo), k = 1,2, by the following formula, 

x k = ip{v)dv + y k 
Jo 

Vk 

(exp(2a(v - v)/a 2 ) - l)dv + y k 



V — Vk 

(-(cr 2 /2a)exp(2a(u - v)/a 2 ) - v) 



+ Vk 

v=0 



= {a 2 /2a){exp(2av /(T 2 ) - exp(2a(t> - y k )/o 2 ))- (3.12) 
Choose yi < j/3 < j/4 < 2/2 and t>i < «o such that, 

af^G Vl (v)dv f£<p(z)dz 



- 1 , ry 2 / \ 1 c o- ( 3 - 13 ) 

2/2 - yi + J yi ip{z)dz V2-y\+ J yi ip{z)dz 

Recall that [a] denotes the smallest integer greater than or equal to a. Let [a\ 
denote the largest integer smaller than or equal to a. 

Let ci = 1 — pi and P2 = 1 — C1/8. By the law of large numbers, we can find a 
large to and make £0 > smaller, if necessary, such that if t > to and e G (0, eo) 
then with probability greater than P2, the number of processes in the interval 
[0, yi] is smaller than or equal to (a/e) J Q V3 G Vl (v)dv. If this event holds then there 
are exactly \(a/e) J^ 3 G Vl (v)dv \ processes Zq in some (random) interval [0,2/5] with 
2/5 > Vx- This implies that there are exactly [(a/e) J^ 3 G vi (v)dv~\ processes Xq in 
[0, £ [(a/e) Jq 3 G Vl (v)dv~\ + 2/5]. For fixed ?/i and j/3, we make V\ < vo larger, if 
necessary, so that 



■V3 

(a/e) / G Vl {v)dv 



rV3 rV3 
+ 2/5 > a / G Vl (v)dv + y 5 > a / G Vl (v)dv + y x 
Jo Jo 
rVi rVi 
>a G Vo {v)dv + yi= / <p(v)dv + yi = x\. 
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Hence, the number of Xq's in the interval [0,.ti] is smaller than or equal to 



(a/e) J Q V3 G Vl (v)dv, with probability greater than p 2 . 

We can make to larger and £o > smaller, if necessary, so that by the law of large 
numbers, if t > io and £ £ (0, £o) then with probability greater than p 2l the number 
of processes Z^ s in the interval [0, 2/2] is greater than or equal to (a/e) J Q V4 G Vl (y)dv. 
If this event holds then there are exactly \_(a/e) J Q y4 G Vl (v)dv\ processes Zq in 
some (random) interval [0,2/6] with 2/6 < 2/2- This implies that there are exactly 
[(a/e) J a V4 G Vl (v)dv\ processes Xq in [0,£ \(a/e) J^ 4 G Vl (v)dv\ + y^}. Note that, 



+ 2/6 < a / G Vl (v)dv + y 2 < a / G Vo (v)dv + y 2 = x 2 . 



rV4 

(a/e) / G Vl (v)dv 
Jo 

Hence, the number of Xq's in the interval [0, X2] is greater than or equal to 
(a/e) J Q V4 G Vl (v)dv, with probability greater than p 2 . 

Let d be defined as in (3.9) but relative to / in place of I k . The two events 
described in the last two paragraphs hold simultaneously with probability greater 
than 1 — ci/4. Then the number of Xq's in [s^aia] is greater than or equal to 
(a/e) l v y l G Vl (v)dv. This and (3.13) imply that 



d t ([xi,x 2 }) > 



e [(a/e) G Vl (v)dv) _ a ^ G Vl (v)dv 



x 2 - xi J Q V2 ip(z)dz + 2/2 - J Q V1 (p(z)dz - 2/1 



2/2 - 2/1 + J V1 <p{z)dz 2/2 - 2/i + J yi <p(z)dz 
It follows from (3.12) that 

cxp(2a(vo — 2/2 )/°' 2 ) = cxp(2av /a 2 ) — 2ax 2 /a 2 
and, therefore, for v < y 2l 

<fi(v) = exp(2a(«o — V )/ <T ' 2 ) — 1 > exp(2a(fo — y 2 )/cr 2 ) — 1 

= cxp(2awo/cr 2 ) — 2ax 2 /a 2 — 1. 

We combine this estimate with (3.14) to see that, with probability greater than 
l-ci/4, 

d t ([x 1 ,x 2 ]) > , r y 2 . r: - c 

2/2 - 2/1 + J yi (p(v)dv 

/^ 2 (exp(2aw /o- 2 ) - 2aa;2/cr 2 - l)dv 
~ 2/2 - 2/1 + /y i 2 (exp(2au /cr 2 ) - 2ax 2 /a 2 - l)dv 
_ (2/2 - 2/i)(exp(2at;o/a 2 ) - 2ax 2 /ff 2 - 1) 
2/2 - 2/1 + (2/2 - 2/i)(exp(2aw /cr 2 ) - 2ax 2 /a 2 - 1) 
exp(2ai>o/a 2 ) — 2ax 2 /a 2 — 1 
exp(2awo/c 2 ) — 2ax 2 /a 2 

(cr 2 /2a)(exp(2aw /a 2 ) - 1) - x 2 
~ (a 2 /2a)(cxp(2av /(T 2 ) - 1) - x 2 + a 2 /2a 

= i ,^2/0 -co- (3.15) 

1 — x 2 + a z /2a 
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The last equality follows from (3.11). 

Recall that n t is the number of processes Z k alive at time t. Note that for any 
< t\ < t<i < oo with t2 — t\ > e/a, we have, 

n t2 ~ THi < (o/e)(*2 - (3-16) 
Fix arbitrary t\ > to and choose 5 > such that 



G Vl (v)dv + 5< / G VQ (v)dv. (3.17) 
Jo 

We make £o > smaller, if necessary, so that, by the law of large numbers, if 
e G (0, £o) then with probability greater than 1 — ci/4, for all of the form 
s k = kS/2, k = 0, . . . , \t 1 /5\ + 1, we have n Sk < (a/e) (j^ 1 G Vl (v)dv + 5/2). It 
follows from (3.16) and (3.17) that for e < aS/2, 



sup n t < (a/e) I / G Vl (v)dv + 5 \ < (a/e) / G VQ (v)dv. 
0<t<t! \Jo J Jq 

Suppose that this event holds. Then, for every t < ti, the right edge of the rightmost 
interval 7 t is to the left of 



e(a/e) / G Va (v)dv +Vi = G Vo (v)dv + v + (vi - v ) = 1 + («x - u ) < 1. 
Jo Jo 

(3.18) 

The second equality in the above formula follows from (3.11). 

Recall the definitions given before the statement of the theorem. A process Z k 
is killed when the right end of the rightmost interval I k hits 1. Since the processes 
Z k are driven by the same Brownian motions as Z k , (3.18) implies that every Z k 
has a longer lifetime than Z k . This implies that ([xx, x%\) > d tl ([xi,X2])- We 
combine this with (3.15) to conclude that, with probability greater than 1 — Ci/2, 

1—X2 

l-x 2 + a 2 /2a 

A completely analogous argument shows that, with probability greater than 1— c\/2, 

d tl ([2:1,0:2]) < z Xl 2/ +c . 

1 — x\ + cr /2a 

This completes the proof of the theorem. □ 



d tl ([xi,x 2 ]) > d tl ([xi,a;2]) > — ; — jToT ~ c o- 



4. Discussion 

Remark 4.1. In the following remarks we will refer to the model analyzed in Sec- 
tion 3 as model (C). We will present another model, which we will call (R). Here, 
C represents the "constant" rate of influx of new particles, and R stands for the 
"random" rate of influx. Model (R) consists of a constant number n of "particles" 
I k which are confined to the interval [0,1]. The fc-th leftmost "particle" is rep- 
resented by an interval I k = (X k - e/2,X k + e/2). The intervals I k and P are 
always disjoint. The processes X k are driven by independent Brownian motions 
with the diffusion coefficient a 2 . When two intervals I k and P collide, they reflect 
instantaneously. The number of particles n is such that ne = b, a constant. When 
X k hits 1, it jumps to 0. We conjecture that as e — > 0, the mass density d in the 
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stationary regime for this process has the density (1— x)/(l — x + a 2 /(2a)), just 
like in model (C), where a, a and b are related by the following formula, 



Heuristically, we expect processes X k in model (R) to jump at a more or less 
constant rate in the stationary regime, so this is why we believe that models (R) 
and (C) have the same hydrodynamic limit. We chose not to analyze model (R) 
in this paper as it appears to be harder from the technical point of view while it 
seems to illustrate the same phenomenon as model (C). 

Remark 4.2. Model (R) is closely related to a model studied by T. Bodineau, 
B. Derrida and J. Lcbowitz (Bodineau et al., 2010). In their model, one considers a 
periodic system of L sites with N particles. The particles perform random walks but 
cannot cross each other — it is the symmetric simple exclusion process. At some fixed 
edge, the jump rates are no longer symmetric but jumps occur with rates p in one 
direction and 1 — p in the other direction. The case p = 1 corresponds to model (R) 
described in the previous remark. In the stationary state, the rescalcd density varies 
linearly on the unit line segment, with a discontinuity located where the jump rates 
are biased. Hence, away from the singularity, the stationary empirical distribution 
is harmonic for the generator of the single particle process, i.e., Laplacian. This 
does not apply to the density of mass d in our models (C) and (R). 

Remark 4.3. Since the density d of the intervals I k has the form (1 — x)/(l — 
x + a 2 /(2a)), it is elementary to check that the typical gap size between I k, s is 
ea 2 /(2a(l — x)). In a model with infinitely small particles X k , the gap size is 
also c/ (1 — x) but we do not have any heuristic explanation why the two functions 
representing the typical gap size should have the same form in both models. 

Remark 4.4. We conjecture that the motion of an individual tagged particle I k 
in model (C) converges, as e — > 0, to a deterministic motion with fluctuations 
having the "fractional Brownian motion" structure. In other words, we conjecture 
that the fluctuations are Gaussian with the local scaling of space and time given by 
Ax = (At) 1 / 4 . Our conjecture is inspired by the results in Diirr et al. (1985); Harris 
(1965); Swanson (2007, 2008) on families of one dimensional Brownian motions 
reflecting from each other. 

Remark 4.5. A d-dimensional counterpart of model (C) can be represented as fol- 
lows. Let I k be balls with radius e and center X k . Our d-dimensional model 
consists of a constant number n of / fe 's which are confined to the cube [0, l] d . The 
balls I k and P are always disjoint. The processes X k are driven by independent d- 
dimcnsional Brownian motions with the diffusion coefficient a 2 . When two balls I k 
and P collide, they reflect instantaneously. Let Si and S r be two opposite (d— 1)- 
dimensional sides on the boundary of [0,l] d . Balls I k are pushed into the cube 
through Si at a constant rate a, i.e., e - ^ -1 ) balls are pushed into the cube every 
e/a units of time, uniformly over Si. Once inside the cube, the balls reflect from all 
sides except S r . When a ball hits S r , it is removed from the cube. We conjecture 
that in the stationary regime, when e is small, the density of the mass analogous 
to d will be a function of the distance x from Si, i.e., a function of depending only 
on one coordinate. We do not see any obvious reason why the density should have 
the form (1 — x)/(l — x + o 2 /(2a)). In relation to Remark 4.4, we conjecture that 
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the motion of a tagged particle in the present model is diffusive, with the diffusion 
coefficient depending on x. If this is true, it means that the "pressure" applied to 
particles in one direction can have a dampening effect on the size of oscillations of 
an individual particle in orthogonal directions. 
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